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Limit on Lorentz and CPT Violation of the Bound Neutron Using a Free Precession 

3 He/ 129 Xe Co-magnetometer 

C. Gemmel, 1 W. Heil/'E S. Karpuk, 1 K. Lenz, 1 Yu. Sobolev, 1 K. Tullney, 1 M. Burghoff, 2 W. Kilian, 2 
S. Knappe-Griineberg, 2 W. Miiller, 2 A. Schnabel, 2 F. Seifert, 2 L. Trahms, 2 and U. Schmidt 3 

1 Institut fur Physik, Johannes Gutenberg-Universitdt, 55099 Mainz, Germany 
2 Physikalisch-Technische-Bundesanstalt (PTE) Berlin, 10587 Berlin, Germany 
3 Physikalisches Institut, Universitdt Heidelberg, 69120 Heidelberg, Germany 
(Dated: January 20, 2011) 

We report on the search for Lorentz- violating sidereal variations of the frequency difference of co- 
located spin species while the Earth and hence the laboratory reference frame rotates with respect 
to a relic background field. The comagnetometer used is based on the detection of freely precessing 
nuclear spins from polarized 3 He and 129 Xe gas samples using SQUIDs as low-noise magnetic flux 
detectors. As result we can determine the limit for the equatorial component of the background 
field interacting with the spin of the bound neutron to be b± < 3.7 • 10 -32 GeV (95% C.L.). 
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A great number of laboratory experiments have been 
designed to detect diminutive violations of Lorentz in- 
variancc. Among others, the Hughes-Drever-like experi- 
ments [HQ have been performed to search for anomalous 
spin coupling to an anisotropy in space using electron 
and nuclear spins with steadily increasing sensitivity 0- 
Lorentz- violating theories should generally predict 
the existence of privileged reference systems. In contrast 
with the situation at the end of the 19th century, we have 
a rather unique choice nowadays for such a "preferred in- 
crtial frame," i.e., the frame where the Cosmic Microwave 
Background (CMB) looks isotropic. Trying to measure 
an anomaleous coupling of spins to a relic background 
field which permeates the Universe and points in a pre- 
ferred direction in spacetime as a sort of new aether wind 
is a modern analogue of the original Michelson-Morlcy 
experiment. 

The theoretical framework presented by Kostelecky 
and colleagues parametrizes the general treatment of 
CPT- and Lorentz violating effects in a standard model 
extension (SME) The SME was conceived to fa- 

cilitate experimental investigations of Lorentz and CPT 
symmetry, given the theoretical motivation for violation 
of these symmetries. Although Lorentz-breaking inter- 
actions are motivated by models such as string theory 
[H, E3]> loop quantum gravity [ilj-EHj]. etc. (i.e., fun- 
damental theories combining the standard model with 
gravity) , the low-energy effective action appearing in the 
SME is independent of the underlying theory. Each term 
in the effective theory involves the expectation of a tensor 
field in the underlying theory. These terms are small due 
to Planck-scale suppression and, in principle, are mea- 
surable in experiments. Predictions for parameters in the 
SME for a loop quantum gravity system with a preferred 
frame were discussed, e.g., in Ref. [22| . 

The SME contains a number of possible terms that 



couple to the spins of standard model particles like the 
electron, proton, and nucleon (mostly the bound neu- 
tron) f23[. These terms have set the most stringent lim- 
its on CPT and Lorentz violations. To determine the 
leading-order effects of a Lorentz violating potential V, 
it suffices to use a non-relativistic description for the par- 
ticles involved given by [23j | 



V = -b 



(with J = X, Y, Z ; w — e,p,n) . (1) 



The most sensitive tests were performed using a 3 He- 
129 Xe Zeeman maser to place an upper limit on the 
neutron coupling to the anomalous field of 6™ = 

~ b x) 2 + (by) 2 < 10 -31 GeV @,[T^ and, recently, by 

use of a K- 3 He co-magnetometer thereby improving the 
previous limit by a factor of 30 [3] . An essential assump- 
tion in these so-called clock comparison experiments is 
that the anomalous field bj does not couple to magnetic 
moments but directly to the sample spins a J . This direct 
coupling allows comagnetometry that uses two different 
spin species to distinguish between a normal magnetic 
field and an anomalous field coupling. 

The comagnetometer used for the presented measure- 
ments is based on the detection of freely spin process- 
ing nuclear spins from polarized 3 He and 129 Xe samples 
gas with SQUIDs as low-noise magnetic flux detectors. 
Like in P, we search for sidereal variations of the 
frequency of colocated spin species while the Earth and 
hence the laboratory reference frame rotates with respect 
to a relic background field. The observable to trace pos- 
sible tiny sidereal frequency modulations is the combina- 
tion of measured Larmor frequencies given by 
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By that measure the Zeeman term is eliminated and 
thus any dependence on fluctuations and drifts of the 
magnetic field. For the 3 He/ 129 Xe gyromagnetic ratios 



2 



we took the literature values [13, HE] given by ^He/lXe = 
2.75408159(20). 

The essential difference, in particular, from [H, [l(|, is 
that by monitoring the free spin precession, an ultrahigh 
sensitivity can be achieved with a clock which is almost 
completely decoupled from the environment. The de- 
sign and operation of the two-species 3 He/ 129 Xe comag- 
netometer has been shown recently [26| . Briefly, in our 
measurements, we used a low-T c DC-SQUID magnetome- 
ter system inside the strongly magnetically shielded room 
BMSR-2 at PTB [27|. A homogeneous guiding magnetic 
field Bq of about 400 nT was provided by one of the two 
square coil pairs which were arranged perpendicular to 
each other in order to manipulate the sample spins, e.g., 
7r/2 spin flip by nonadiabatic switching. The maximum 
field gradients were about 33 pT/cm. The 3 He/ 129 Xe nu- 
clear spins were polarized outside the shielding by means 
of optical pumping. Low-relaxation spherical glass ves- 
sels (R=3 cm) were filled with the polarized 3 He/ 129 Xe 
gases and placed directly below the Dewar as close as pos- 
sible to SQUID sensors, which detect a sinusoidal change 
in magnetic flux due to the spin precession of the gas 
atoms in the glass cell. In order to obtain a high common 
mode rejection ratio, gradiometric sensor arrangements 
are commonly used. For our analysis it was sufficient to 
use a first-order gradiomctcr in order to suppress envi- 
ronmental disturbance fields. 

Nitrogen was added as buffer gas to suppress the van 
der Waals spin relaxation of 129 Xe (28|. In the regime 
of motional narrowing, i.e., at gas pressures of order 
mbar and at low magnetic fields [29|, |3(| , transverse spin- 
relaxation times T 2 * of up to 60 h have been measured 
for 3 He. The actual limitation in the Xe of xenon is 
given by the relatively short wall relaxation time of 8 h < 
Tj^. aU < 16 h. Therefore, the total observation time T of 
free spin-precession of our 3 He/ 129 Xe comagnetometer is 
set by this characteristic time constant. According to the 
Cramer-Rao Lower Bound (CRLB) [3l|, the accuracy by 
which the frequency of a damped sinusoidal signal can 
be determined is given by 

g ^ (2,).5^^.T 3 / 2 X V / ^^- ( 3 ) 

SNR denotes the signal-to-noise ratio, fsw the band- 
width, and C(T, T 2 ) describes the effect of exponential 
damping of the signal amplitude with T 2 . For observa- 
tion times T < T 2 *, C(T, T 2 *) is of order one. Deviations 
from the CRLB power law, due to noise sources inher- 
ent in the comagnetometer itself, did not show up in Al- 
lan standard deviation plots used to identify the power- 
law model for the phase noise spectrum of our runs with 
T « 14 h, typically (2fj|. 

The recorded signal is a superposition of the 3 He and 
129 Xe precession signals at Larmor frequencies whc = 
7Hc • B a ks 2ir ■ 13.4 Hz and w Xc = 7x c • B « 2n ■ 4.9 Hz 
as shown in FIG. [TJ Analogue to similar problems of 
data analysis [32| phases of subdata sets were analyzed: 
The data of each run (j = 1, ... ,7) were divided into 
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FIG. 1. Typical subdata set of 3.2 s length showing 
the recorded SQUID gradiometer signal of the precessing 
3 He/ 129 Xe sample spins (sampling rate: r a =250Hz). The 
uncertainty at each data point is ± 34 fT (k=l) and therefore 
less than the symbol size. The signal amplitudes at the begin- 
ning of each run were typically She ~13 pT and Sx e ~4 pT. 

sequential time intervals (i) of r = 3.2 s (i = 1, . . . ,Nj). 
The number of obtained subdata sets laid between 13350 
< Nj < 18000 corresponding to observation times Tj of 
coherent spin precessions in the range of 12 h < Tj < 16 h. 
For each subdata set a x 2 minimization was performed, 
using the fit function 

A\t) =A\, C ■ sin (c4 e i) + B l Hc ■ cos (w^f) + 
^kc • sin (w Xe i) + B Xo • cos (w Xc t) + 
(4 + c\ iD ■ t) (4) 

with a total of 8 fit parameters. Within the relatively 
short time intervals, the term (cq + cj in • t) represents 
the adequate parameterization of the SQUID gradiome- 
ter offset showing a small linear drift due to the elevated 
1// noise at low frequencies (< 1 Hz) 26]. On the other 
hand, the chosen time intervals are long enough to have 
a sufficient number of data points (800) for the x 2 min- 
imization. The sum of sine and cosine terms are chosen 
to have only linear fitting parameters for the subdata set 
phases which are given by 

ip* = arctan (B^A 1 ) . (5) 

The normalized % 2 (% 2 /d.o.f) of most subdata sets (i) is 
close to 1 which is consistent with the assigned uncer- 
tainty to each data point of ± 34 fT (fc=l); see FIG. Q] 
The latter value is the typical noise signal N s derived 
from the mean system noise p s ~ 3 fT/\/Hz in the 
recorded effective bandwidth of 100 Hz. Jumps in the 
SQUID signal in the order of 1 pT caused by external 
disturbances gave x 2 /d.o.f^>l for the respective sub- 
data sets. In the analysis we therefore disregarded all 
subdata sets with x 2 /d.o.f>2 (< 0.5% in total). For 
each subdata set of chosen time interval -1.6 s< (t — 
h—i,j) < +1-6 s (see FIG. [1]), we finally obtain numbers 
for the respective fit parameters w^ e , w Xo , <^ c , </? Xe in- 
cluding error bars. 

In a further step, we can deduce values for the average 
frequency uji — X)i=i w * f° r each run. The accumu- 
lated phase (omitting the index j) is then determined to 
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beE] 

$ (t = tut) = $ (t = (m - l)r) + w • r + </? m 

- mod [$ (t = (m - 1)t) + cjt; 2tt] (6) 

with m=l,...,N-l and = 0) = ip 1 being the phase 
offset of the first time interval. Following Eq. (f2]) the ex- 
tracted phase difference A$ (1) (t) = $£j (t) - (7h /7x c ) • 

^Xe (*) i s plotted for run f . A$W (t) is expected to 
be constant if there is no sidereal modulation of the 
spin-precession frequency and/or no other drifts and 
noise sources. Nevertheless, in addition to an arbitrary 
phase offset an almost linear time dependence is seen in 
FIG. HJa). The dominant contribution is caused from 
the Earth's rotation, i.e., the rotation of the SQUID de- 
tector with respect to the processing spins. For the loca- 
tion of the PTB Berlin, Germany (0=52.5164° north) and 
the angle between north-south direction and the guiding 
magnetic field (p=28°), the linear term in the weighted 
phase difference due to Earth's rotation is given by [26[ 



$Earth = -f^SD ' (1 - 7Hc/7Xg) ' COS p ■ COS • t 

= 6.87263 x l(T 5 rad/s • t . 



f2sD is the angular frequency of the sidereal day with 
fl SD = 2it/T sd = 27r/(23 h : 56 min : 4.091 s ). Subtract- 
ing this term from A<fcW(i), we get the corrected phase 

A$corr(t) which is plotted in FIG. [2£a), too. Let us as- 
sume that there is no sidereal variation of the 3 He/ 129 Xe 
frequencies induced by Lorentz-violating couplings, and 
then A<I>cOT r (i) can be described best by 



*g (t) = 



A42(*-*o,i) 



-E 



U) 

He 



+E 



Xe 



exp 



exp 




(8) 



with 



for toj <t< (toj 
elsewhere. 



toj is the absolute starting time of each run. Our inter- 

(i) 

pretations of the terms are as follows: $0 is a general 
phase offset and Awj^(t — toj) is an additional linear 
phase shift mainly arising from deviations of the gyro- 
magnetic ratios of 3 He and 129 Xe from their literature 
values due to chemical shifts and uncertainties in the sub- 
traction of $Earth [26| . The two exponential terms with 
(7) (i) 

amplitudes and E^ reflect the respective phase 



1 As the maximal frequency deviation Au from the mean ui J was 
smaller than 5 ■ 10 — 6 rad/s in the course of one run [26(1 . we had 
at all times Acu ■ t <S 2tt. 
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FIG. 2. (a) Measured phase differences A&^ft) for run 1 and 
the corresponding corrected phases A$corr(i) after subtraction 
of the effect of the Earth's rotation, (b) Phase residuals after 
subtraction of phase drifts given by the fit model of Eq. ([8} 
(one data point comprises 20 subdata sets, i.e., At = 64 s). 



shift due to demagnetization fields in a nonideal spherical 

cell seen by the spin ensembles (self-shift). These phase 

(7) 

shifts are directly correlated to the decay times T 2 h and 

x e °f the respective signal amplitude of the precessing 

helium and xenon spins (26|. As the T* {]) times can 
be determinded independently for both spin species from 
the experiment, four fit parameters are left for each run, 
such that the fit model is basically a linear function in 
parameters. 

Fitting the corrected phase difference A$i J o„(t) to 
Eq. (|SJ) and subtracting the fit function from A$id rr (t) re- 
sults in the phase residual as shown for run 1 in FIG.[U(b). 
Because of the exponential decay of the signal ampli- 
tudes, mainly that of xenon with the shorter T^Xc °f 
only 4-5 h, the SNR decreases resulting in an increase of 

the residual phase noise, i.e., cr$ ircs cx exp ^/^2,Xe) [111- 

In the last step, a piecewise fit function was defined, 
which is a combined fit to all seven runs, now including 
the parameterization of the sidereal phase modulation 



4 




140 



FIG. 3. Corrected phase differences A<&corr(t) with combined 
fit function c&f] 3 ^) (white solid line) for the seven runs (one 
data point comprises 20 sub-data sets, i.e., At = 64 s). In 
order to present these results in a common plot, the general 
phase offset was subtracted from A$io TT (t) for each run. 



(*)=5>g>(*) 

+{ a s ■ sin (fi S D it - to,i) + </?sd) 

- a c ■ cos (fi S D (t - t 0>1 ) + ipso) } (9) 



ipSD = 27r-isD represents the phase offset of sidereal mod- 
ulation at the local sidereal time (vernal equinox J2000.0) 
on 2009 March 21 at 20:52 UT (universal time) which is 
the starting time io.i of the first run. Neglecting multi- 
ples of 24 h the local sidereal time is 9.7 h, which in units 
of sidereal day gives isD = 0.4053. 

Figure |3] shows the corrected phase differences $corr(X) 
together with the fit function (t) (white solid line) for 
all seven runs. The x 2, /d.o.i of the fit gave 1.868, which 
shows that the phase model of Eq. © may be somewhat 
incomplete or, what is more likely, the phase errors are 
underestimated in the analysis of the subdata sets. In 
order to take an (unknown) uncertainty into account, the 
errors on the phases were scaled to obtain a x 2 /d.o.f of 
one, as recommended, e.g., by Refs. [33|,|34j]. In Table „ 
(2nd row) the fit results for the amplitudes a c and a s of 
the sidereal phase modulation are shown together with 
their correlated and uncorrelated la errors. 

It is noticeable that the uncorrelated error which rep- 
resents the integrated measurement sensitivity of our 
3 He/ 129 Xe comagnetometer is about a factor of 50 less 
than the correlated one. The big correlated error on a s 
and a c is caused by a piecewise similar time structure of 
$1+ (t) and the sidereal phase modulation in the fit func- 
tion of Eq. ©. On a closer look, this can be traced back 

to the relatively short T 2 *xe times (compared to Tsd) that 
enter in the argument of the exponential terms of Eq. © . 
Therefore, the present sensitivity limit of our 3 He/ 129 Xe 
comagnetometer is set by the correlated error. In order to 
substantiate that more clearly, we changed the fit model 
of Eq. © by taking multiples of fisD (^sd = <7"^sd), i-e., 
replacing Tsd by T SD = Tsd / <?■ The results show that the 
correlated error approaches the uncorrelated one already 



for g > 3 (see Table^. The uncorrelated error, however, is 
only marginally affected by this procedure, as expected. 
From Table n (2nd row) we now extract the rms value of 
the sidereal phase amplitude $sd = -\/a 2 + a 2 , yielding 
(2.25 ± 2.29) mrad (95% C.L.). This result is consistent 
with the abscence of Lorentz- and CPT- violating effects, 
giving reasonable assumptions about the probability dis- 
tribution for $sd [35j . 

In terms of the SME [23| we can express the sidereal 
phase amplitudes according to 



2tt 



9. 



Svx and a c 



SD 



2tt 

^SD 



Suy 



(10) 



with 



2tt|^x,y| -h= 2-(l- 7 Hc/7Xc)-sinx-fe x ,Y ■ (H) 

X is the angle between the Earth's rotation axis and the 
quantization axis of the spins with \ = arccos (cos# • 
cosp) = 57°. Within the Schmidt model (3(|, the valence 
neutron of 3 He and 129 Xe determines the spin and the 
magnetic moment of the nucleus. Thus, our 3 He/ 129 Xe 
comagnetometer is sensitive to the bound neutron pa- 
rameters &xy- F rorn that, we can deduce numbers for 

in 

°X,Y- 

6£ = (3.36 ± 1.72) ■ 10" 32 GeV (la), (12) 
6 Y = (1.43 ± 1.33) ■ 10~ 32 GeV (la), (13) 



< 3.7 • 10" 32 GeV at 



which can be interpreted as 

95% confidence level for the upper limit of the equatorial 
component of the background tensor field interacting 
with the spin of the bound neutron. For the calculation 
of the upper limit on b\ Eq. (fTUf and Eq. (fTTj) were used 
putting in the 95% C.L. for the rms value of the sidereal 
phase amplitude $sd- 

Further improvements for Lorentz and CPT tests 
using the free spin-precession 3 He/ 129 Xe comagnetome- 
ter can be achieved via two mayor steps: First, the 
relatively short wall relaxation time of 129 Xe limiting 
the total observation time T of free spin precession 
has to be increased considerably (T 1:Wa n Tsd) such 
that we approach the measurement sensitivity given by 
the uncorrelated error. Since the latter one follows the 
ocT~ 3 / 2 power law according to CRLB of Eq. d3]), the 
longer observation time T will lead to an additional 
increase in sensitivity. Second, the number of measure- 
ment runs has to be increased to a period of 100 days. 
Besides gain in statistics, the long time span provides 
an important separation between sidereal and possible 
diurnal variations. 

This work was supported by the Deutsche Forschungs- 
gemeinschaft (DFG) under Contract No. BA 3605/1-1. 
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TABLE I. Results for the sidereal phase amplitudes a c and a s together with their 
correlated and uncorrelated la errors (2nd row) determined by a x 2 minimization 
using the fit model of Eq. ([9]). In order to demonstrate the strong dependence of the 
correlated error on the angular frequency of the sidereal day £Isd, corresponding fit 
results are shown for multiples of QsW- ^sd — 9 ' ^bd- 

q c (mrad) a c °™ (mrad) qg c ncorr (mrad) a s (mrad) <? rr (mrad) < lcorr (mrad) 



0.5-fisD 


3.353 


6.572 


0.018 


0.488 


7.991 


0.016 


^SD 


-0.882 


0.814 


0.015 


-2.067 


1.057 


0.019 


2-f2sD 


-0.048 


0.120 


0.016 


-0.149 


0.112 


0.017 




-0.184 


0.052 


0.019 


-0.011 


0.043 


0.016 


4-f2 SD 


-0.001 


0.034 


0.018 


0.057 


0.030 


0.016 
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